Dislocations and Bragg glasses in two dimensions 
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We discuss the question of the generation of topological defects (dislocations) by quenched 
disorder in two dimensional periodic systems. In a previous study [Phys. Rev. B 52 1242 (1995)] 
we found that, contrarily to d = 3, unpaired dislocations appear in d = 2 above a length scale 
£d, which we estimated. We extend this description to include effects of freezing and pinning of 
dislocations at low temperature. The resulting £d at low temperature is found to be larger than our 
previous estimate, which is recovered above a characteristic temperature. The dependence of £d in 
the bare core energy of dislocations is a stretched exponential. We stress that for all temperatures 
below melting £d becomes arbitrarily large at weak disorder compared to the translational order 
length R a 3> a. Thus there is a wide region of length scales, temperature and disorder where Bragg 
glass like behavior should be observable. 
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An outstanding current problem in condensed matter 
physics is to understand the physical properties of pe- 
riodic systems in presence of point impurities. This is 
important for numerous experimental systems, such as 
vortices in type II ^uyperconductorsEm, Wigner crystal- 
lization of electrons3~Q, charge density wavesQ etc.. We 
have proposed^ that, contrarily to previous claims, such 
systems generically possess at weak disorder in three 
dimensions a distinct thermodynamic glass phase, the 
Bragg glass, with perfect topological order (dislocation 
free) and quasi long range translational order. An im- 
mediate consequence of this theory is that a phase tran- 
sition away from the Bragg glass at which dislocations 
suddenly proliferate must occur upon increase of disor- 
der (or field) which is determinant for the phase diagram 
of type II superconductorsa'EI. Thus, whether dislocations 
are generated or not by disorder has emerged as a cen- 
tral question for these systems. In three dissensions there 
is increasing theore±kalt9il3 , numericalll3ila and ex- 
perimental evidenceEjtLj that this topologically ordered 
Bragg glass phase exists. By contrast, the strong disor- 
der glassy state in three dimensional superconductors is 
still poorly known beyond the fact that it must contain 
topological defects. It is yet unclear whether this state 
is a distinct thermodynamic phase or simply a crossover 
from a pinn ed , li quid. The initial proposal of a "vortex 
glass" phaseEjEy (claimed to contain topological defects) 
now clearly cannot stand at weak disorder. It is even un- 
cleaip-|in view of recent numericalE3E3 and experimental 
datacil whether such a phase, as described in Ref. ^0 21 



exists at all in real (i.e with screening) three dimensional 
superconductors. 

In two dimensions it is easier to generate dislocations 
and we have obtainedEl that within the conventional per- 
turbative RG analysis unpaired dislocations appear be- 
yond a length scale £u which we have estimated. A 



e specific case of 
'' Recent numeri- 



similar result was also obtained fc 
two dimensional triangular latticel 

cal studiesHEZl in d = 2 support these predictions. Note 
however that the question of the stability of the Bragg 
glass can only be decided numerically by investigating 
the weak disorder region. Thp RG flow usedcl, although 
based on perturbation theoryEScj, already captured the 
delicate balance between the elastic energy cost and the 
energy gain due to disorder. One striking consequence 
was that although unpaired dislocations are generated, 
the scale at which they appear, obtained as: 



£d ~ R a e 



c(ln(fl Q /a)) 1/2 



(1) 



can be made, at weak disorder R a 3> a, arbitrarily large 
compared to the length R a at which the displacements 
become of the order of the lattice spacing. 




(a) " ' (b) 

FIG. 1. (i) Conventional picture of a solid broken up 
by disorder in domains of size R a with unpaired dislocations 
(black dots) appearing at the same scale, (ii) correct picture 
for weak disorder: the scale £d at which unpaired dislocations 
appear is larger than the scale R a at which translational order 
starts decaying slowly. 
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Below melting we thus obtained that even in rf = 2 
the naive picture of the system as breaking up in crys- 
tallites of size R a is still incorrect as illustrated in Fig. 
ffl. Even if d = 2 is the lower critical dimension of the 
Bragg glass, since can be very large in practice a wide 
regime of effective Bragg glass behaviour should still be 
observable. As pointed out in Rcf. |25| where a theory 
predicting £d{T) near melting was obtained, this may al- 
low to understand why experiments on two dimensional 
superconductorsEj show a sharp change of behaviour with 
characteristics of melting. 

In this paper, we further examine the two dimensional 
problem. Indeed there has been recent evidence that 
at low temperature the conventional RG equations can 
overesi*wftte^h*|-4mportance of dislocations in related 
modelsEj~E3'Ea'tZlc3. Thus our initial estimate ([[]) of £d , 
which should be accurate at intermediate temperatures, 
may become only a lower bound on at low temper- 
ature. We give here a more accurate description of the 
low temperature properties. This is of importance for 
applications to experiments or simulations on classical 
systems, often performed at low temperature but also for 
ground state properties of quantum disordered .systems, 
such as the two dimensional Wigner crystalocj. 

In order to study dislocation problems in d = 2 let us 
restrict for simplicity to the single component scalar XY 
model in a random field, defined in the continuum by the 
partition sum Z = J D<f>e~ H t T and hamiltonian: 



// = / d 2 X —{V4>{X) -T](X)) 2 

2tt 

-£i (x) cos <t>(x) — £2 {x) sin 4>{x) 



(2) 



with r]i(x)r]j(x') = iraSijS(x — x') and Ci(x)Cj( x> ) = 
gSijS(x ~ x') are Gaussian white noises. Configurations 
with vortices (of integer charge qi at x{) are described by 
decomposing <p{x) = <f>sw{x) + 4>v(x) where 4>sw{x) is 
the smooth (spin wave) field and <j>v(x) the vortex con- 
tribution V x V0y(x) = '^ /i qi2TrS(x — Xi). This model 
contains essential ingredients of a variety of elastic dis- 
ordered systemsEl in d = 2. For instance, in the case of a 
lattice, <p = 2iTu/a can be thought of as a displacement 
field u in units of lattice spacing a. Point like impu- 
rities produce a random potential V{x) which couples 
to the density and leads to (||) with g ~ A^/o 4 pro- 
portional to the amplitude of the disorder with Fourier 
component close to Kq = 2-k/cl (called the pinning disor- 
der) and a ~ Ao/47r J 2 a 2 proportional to the long wave- 
length disorder, where A q = V q V- q . Vortices in the 
field </> correspond to dislocations in the lattice and thus 
will generically be called "dislocations" in the following. 
These defects are characterized by a fugacity y, which 
in the bare model is related to the defect core energy 
y = er Ec -l T . Note that the long wavelength disorder a 
is always generated by coarse graining. The bare model 
(Q) corresponds, for the lattice problem, to the simpler 
case where the correlation length of the disorder r/ is of 
the order of the lattice spacing a. Thus in the model 



specifically studied in this paper, the translational order 
correlation lengths R a - such that relative displacements 
u(R a ) — u(0X-r^ a, is of the same order than the Larkin 
OvchinnikoveL] pinning length R c - such that relative dis- 
placements u(R c ) — u(0) ~ rf, and we will implicitly 
equate them in the following. The situation rj < a (thus 
R c <C R a ) will be briefly mentioned at the end. 

Let us summarize the RG analysis which led to the es- 
timate in Ref. |^ of the scale £d beyond which unpaired 
dislocations appear in d = 2 at weak disorder. The RG 
equations for the fugacity of dislocations, the disorder 
and the stiffness were derived by Cardy and Ostlund 
(CO)Eaa. The fugacity of the vortices y satisfies to low- 
est order in y: 



dy 
dl 



(2 



J 
T 



aJ 2 

J^2 



) V 



(3) 



with T m — J/2 is the pure system Kosterlitz-Thouless 
transition temperature and / = InL is the logarithmic 
scale. The pinning disorder renormalizcs as: 



(4) 



4J = 8T„ 



If dislocations are ex- 



up to 0(g 3 ) terms, with r = I 
and B a nonuniversal constant, 
eluded by hand (setting y — 0) there is a transition at 
T = T g between a high temperature phase (T > T g ) 
where the disorder is irrelevant and a low temperature 
glass phase (T < T g ). For T < T g there is a line of 
fixed points g = g* which describes a 2d Bragg glass 
phase where displacements grow as u ~ In x and beyond 
which translational order slowly decays. This asymp- 
totic behaviour is reached at the length R a ~ 0(^)27 
(where go is the bare value) which for weak disorder is 
R a ^> a. The global structure of the RG (together with 
functional RG extensions) suggests that the fixed line 
continues down to T = and that correlations behave as 
{{(j)(x) -|^0)) 2 ) = C{T)\n{x/R a ) with a T dependent 
prefactoiO. It is perturbatively controlled only near 2^-. 
where g* is small, and the correct universal prefactorEJ 
C(T) — 2t 2 was derivedE3 to lpy es t order in Tj—Numer- 
ical sigQ^3ns have measuredcl C(T) near TgC5 and at 

When dislocations are allowed in the model, one finds 
that they are perturbatively relevant above T m leading 
to a liquid. Below T m however, since the bare value of 
a is small for weak disorder, one could naively conclude 
from (||) that dislocations are suppressed at low enough 
temperature. As was discussed in Ref. || this is not the 
case and in fact the 2d Bragg glass fixed point is unstable 
to dislocations. Such an instability is a peculiar feature 
of d = 2 and does not occur in rf = 3 where the Bragg 
glass phase is stable with respect to dislocations at weak 
disorder. The instability occurs because in rf = 2 the 
long wavelength disorder a is also renormalized: 



& a A 2 



(5) 
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Tg) and thus 



- a' + 

2 _ a_2 



to lowest order (with A = B 2 at T 
grows unboundedly with the scale as a(l) 
2C(T)]n(L/R a ) for L > R a , with 2C(T) = Ag*' = 4r 
near T g . As can be seen from (|J) y(l) starts increasing 
beyond a certain length scale and eventually becomes of 
order y{lo) ~ 1 at a scale Id = ln(£rj/a). At that scale 
unpaired dislocations dominate the behaviour and trans- 
lational order is exponentially destroyed beyond £d- 




FIG. 2. Dependence of the dislocation fugacity y in the 
length scale for weak disorder. It is first strongly renormalized 
downwards before it eventually shoots up again and reaches 
values of order unity at £d S> R a - 

However, £d can be very large and in particular much 
larger than R a . This is because for weak disorder and 
T < T m the fugacity of dislocations is first strongly renor- 
malized downwards for scales smaller than R a (see (||)), 
as showijtiin Figure [| Integrating the flow led us to the 
estimated: 



-l)\n(R a /a) 



(6) 



where T m (a) is the boundary of the XY phase in the 
absence oi-pinning disorder (g = 0) as given by the CO 
equationsE3 (see Fig. | below). The case of triangular 
lattices in presence of disorder was studied in more de- 
tails recently using a N = 2 component model, and a 
generalization of both the above CO equations (||4|j|) 
and the KTNHY equations r^wtich describe the fusion of 
pure crystals) was obtainedEacS. It confirms that a simi- 
lar estimate as (|J) holds for lattices and yields a comple- 
mentary formula for £d around the pure crystal melting 
transition. 

To obtain (0) the only assumption made in Rcf. || and 
supported by the structure of the RG flow, was that be- 
yond the scale R a pinning disorder g has reached a fixed 
point and that o~(l) grows as I = bxL. This is equivalent 
to use beyond R a the effective "random stress" model 
which reads in the continuum: 

(fx^-Wix)-^)) 2 (7) 

to which one must add vortices as discussed above. The 
disorder r]i(q)r)j(-q) = with er(<?) ~ ln(l/g), 



H 



now depends on the scale in a logarithmic way as o~{l) ~ 
a* + 2C(T) bx(L/R a ), where C{T) is th^amplitude de- 
fined above. A recent numerical workEZI comes as ad- 
ditional support that this assumption is indeed valid. 
When a (I) = a is scale independent^ this model reduces 
to the "random phase shift model"t3. The advantage of 
model (0) is that it can be treated by RG or by qualita- 
tive arguments even in presence of dislocations. 

Let us now reexamine ^ at low temperature. Be- 
cause of the reentrance of the disordered phase present 
in the CO equationsE3 an extrapolation of (g) would lead 
to a small £d at low temperature (below T~(er) in Fig. 
|]). However, the original CO equations tend to over- 
estimate the effect of dislocations as they, neglect non 
thermalization and pinning of dislocationsa. Indeed, re- 
cent reexamination of the phase diagram of the random 
phase shift ija^del (@) has shown that the reentrant disor- 
dered phaseEHES disappears when these effects are taken 



into accountO. Although in the high temperature region 
T > T m /2 the RG equation (g) was found to be cor- 
rect, new physics arises below the line T < T*(a). In 
order to reexamine our previous estimate (|^) for £q at 
low temperature, we extend to the present case (where 
a (I) depends logarithmically on L) the modified RG 
analysis pcecently developed for the random phase shift 
modeSErMi 

To study the relevance of dislocations it is first useful 
to consider a single dislocation or dipole at T = 0. The 
simplest energy argument, presented first in Ref. [?1 49, 
estimates when it is favorable to place one vortex or a 
dipole in the system of size L. This vortex sees a 2d ran- 
dom potential V(x) with logarithmic correlations, and 
one must thus estimate the minimum energy E m i n of 
this random potential. A reasonable approximation is 
obtained by neglecting correlations but keeping the cor- 
rect local variance aafj: 



1 

L 2 " 



dV 



\fA-KaJ 2 In L 



exp(- 



V 2 



4a J 2 InL 



(8) 



This estimate leads to E m i n ~ —^/8aJlnL. The more 
accurate methods vdaich are now available to take corre- 
lations into accounioEj confirm that the prefactor is ex- 
act, and estimate the (large) corrections to scaling. The 
energy gain by disorder thus overcomes the elastic cost 
E e i = J InL for a > a c = 1/8. Below this value an 
XY phase, dislocation free at. large scale, exists contrar- 
ily to the earlier conclusiorO based on (||). Thus these 
low temperature effects strongly reduce the relevance of 
dislocations. In the random phase shift model, the topo- 
logically ordered Bragg glass phase is thermodynamically 
stable at low temperatures. 

One can thus expect a similar reduction of the impor- 
tance of dislocations in the presence of pinning disorder, 
i.e for the random stress model with a(l) ~ InL com- 
pared to our previous estimates based on the CO equa- 
tions. A straightforward modification of the above ar- 
gument (0), replacing a by a (I) in (g) leads to E min ~ 
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y'Scrp) J In L ~ (InL) 3 / 2 . The energy gain due to disor- 
der now always overcomes the elastic cost J In L at large 
scale and dislocations are always generated in model (0). 
This modification of the argument of Ref. thus allows, 
to recover in a very simple way the (lnL)^ estimatecll 
for the optimal energy of a dislocation pair. 

Howescer these types of energy arguments or the 
analysiscZI for a single dislocation (or dipole) does not 
by itself allow to compute the length scale £d at which 
vortices (dislocations) destroy XY (positional) order. In- 
deed, to destroy the order exponentially one needs a finite 
density of defects. In particular it would be incorrect to 
identify £d as the length scale at which the disorder en- 
ergy becomes of the order of the elastic energy. This 
is clear for instance when looking at the KT transition 
were elastic energy is J In L and entropy is 2T In L. When 
T > T m — J/2 balancing only these two terms would in- 
correctly predict unpaired dislocations of size a near T m , 
when in reality they occur at a much larger length scale 
£ m ~ ae _c ( T-Tm ) 1/2 . To get the correct result one must 
take into account both the defect core energy E c and the 
screening by smaller dipoles. The same effects must be 
also taken into account for the random stress model in 
order to quantify the importance of dislocations. 

We now estimate the length scale £d at low temper- 
ature using a RG analysis for the random stress model 
(Q). In a recent workup it was shown that to describe 
the correct physics in the model with constant a(l) = a, 
one must follow the full distribution of dislocation core 
energies Pi(E c ), found to satisfy a non linear RG equa- 
tion. At T = 0, this non linear RG equation solves the 
problem of energy minimization iteratively on successive 
logarithmic scales. Its solution Pi(E c ), develops broad 
tails at low temperature, as illustrated in Fig. |[ Very 
schematically, while the center of the distribution is bi- 
ased towards the right as ~ JlnL, the width grows as 
~ aJ 2 lnL. For large enough a (a > 1/8) the width 
"wins" and the probability of finding a site with negative 
effective core energy eventually increases with the scale. 
This analysis can be extended in an accurate manner to 
the present case, by studying the non linear RG equation 
in presence of a scale dependent a (I) ~ I. 

For our present purpose it is enough to use the fol- 
lowing simplified version where Pi(E c ) at scale I is given 
by: 



Pl(E c ) ~exp(2Z- 



(E c -Jl- E° c ) 2 
4J 2 /' a(l')dl' + c 



(9) 



This amounts to study the. solution of the linearized ver- 
sion of the RG equationEa which should be a good ap- 
proximation in the small E c far tail of the distribution 
of Pi(E c ) needed here. The distance between unpaired 
dislocations corresponds to the length scale I = Id = 
ln(£o/a) at which Pi(0) ~ 1. Physically this corresponds, 
in the renormalized model, to putting one dislocation per 
site. The above form, for constant a, clearly yields a dis- 
location free phase (Id ~ +oo) for a < 1/8 and prolifer- 



ation of dislocations for a > 1/8. In that case differen- 
tiating Pi (0) from (H) yields back the RG equatiorH for 
the "effective dipole fugacity" y 2 as: 



dy / 1 s 



T <T* = 2a J 



(10) 



identified as y 2 ~ Pi(0) 2 below the line T = T* = 2a J 
where the freezing does occur (see Figure [|). 



P ( E c) 




typ 

FIG. 3. scale dependent distribution Pi(E c ) of the effec- 
tive core energy as discussed in Ref. [}6| The bulk of the dis- 
tribution is centered around a typical value E typ (l) but broad 
tails develop at low temperature 

In the present case the above equation (|]) is comple- 
mented by the equation for a(l) which arises as the solu- 
tion of (U,||). Denoting R a = ae la , and using the simpli- 
fied form g(l) = g* e 2T( - l ~ la ^ for / < l a and g(l) = g* for 
I > l a , we obtain: 



a (I) = a + 



Ag*< 
4t 



a(l)^a* +Ag* 2 (l-l a ) I > l a 



(11) 



with a* = ^j— 

4r 



This is valid near T g , but can be ex- 



tended everywhere by replacing Ag* 2 by 2C(T) in (HI). 
Also, in the range of length scales needed here, and far 
from the critical crossover region ao — 1/8, we can ne- 
glect the renormaiization of a and J by dislocations. Us- 
ing (|nj) , (g) leads to: 



£d ~ R a e 



(12) 



with c w 1/ J C(T). This expression holds in the low 
temperature region T < T* = 2a J and for weak disorder 
R a ^> a. Interestingly, it does have a form very simi- 
lar to (^) which is valid in the high temperature region 
T* <T < T m (a) in Fig. |, but for the fact that now the 
disorder ao plays a role analogous to temperature. The 
expression (n3) smoothly connects with our previous es- 
timate (g) at higher temperature. It would be interesting 
to check the predictions (f||l2|) in numerical simulations. 
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A previous numerical wor on the random XY model 
close to T m is indeed consistent with £]j larger than R a 
although no attempts were made to compare the results 
with (|6|). Finally, very near oq = 1/8 we expect that 
will take a critical crossover scaling form as is the casec3 
near T — T m . 




FIG. 4. Different regions of the temperature T, disorder 
a phase diagram. The dashed lines represent true XY phase 
boundaries (T m (cr) and a — a c — 1/8) for the model without 
pinning disorder g = 0. The dotted lines are the old (incor- 
rect) XY boundary T^(er) obtained via CO-RG, and the new 
freezing line T*(u). In presence of pinning disorder g > 0, a(l) 
grows with the scale and eventually dislocations appear. The 
arrowed lines represent the RG trajectories in the two cases 
(low and intermediate temperature) discussed in the text 

Note that it would be tempting to replace directly a 
by a(l) in ( |To| ) and integrate the flow. This is incorrect, 
as can be seen by writing the RG equation for P;(0) us- 
ing (^) . This is a manifestation of the fact that one must 
here follow a distribution which cannot, at low tempera- 
ture, be parametrized by a single fugacity-like parameter 
V- 

In numerical simulations it may be easier to study in- 
termediate disorder but change artificially the bare core 
energy Keeping the core energy dependence in the 
previous analysis of (^J) leads to the following equation 
for l D = ]n(£ D /a): 



0\2 



= 21 D 



(JId + E°) 
U 2 C{T)l 2 D 



Let us close with an overall perspective on the situ- 
ation in d = 2. The above analysis of the flow of pin- 
ning disorder can be improved, e.g. allowing for more 
general RG flow structure or using functional RG tech- 
niques appropriate to zero temperature (where the model 
(||) develops higher harmonics and nonanalyticity). This, 
we expect should not change the results obtained here, 
up to numerical prefactors, as long as the approximation 
introduced in Ref. || by the random stress model holds 
at large scale. This model allows to treat some of the 
conventionally non perturbative effects related to dislo- 
cation freezing by long wavelength disorder. Proving its 
validity everywhere, or going beyond it, is difficult an- 
alytically, and challenging numerically because of large 
corrections to scaling. Thus it still cannot be excluded 
that more surprises may be in store when further non 
perturbative effects (e.g. the effect of pinning disorder g 
on dislocations) are taken into account. The case rf<^ia 
can also be studied. An intermediate "random manifold 
regime" then exists between R c < L < R a ~ R c (a/rf) a 
with a ~ 3. A (functional) RG analysis can be sketched. 
The RG equation for a becomes dio ~ J^k K 2 9k to low- 
est order, while the gx reach a crossover functional form 
beyond L = R c and another, asymptotic one, beyond 
L = R a . Although a{l) starts growing beyond R c , we 
find that it remains small until R a . This indicates that 
the above results should still hold. 

To conclude we have analyzed the question of the gen- 
eration of dislocations by disorder in d = 2. In an 
ealier studyH based on a random stress model approx- 
imation of the Cardy Ostlund equations we had found 
that contrarily to d — 3, in d — 2 unpaired disloca- 
tions appear beyond a length scale £d. Recent numer- 
ical simulation d 1 ^! 27 ! have reached a similar conclusion. 
At low temperature we have improved our estimate for 
£d taking into account freezing of dislocations by long 
wavelength disorder. Taking into account these effects 
increases compared to our previous estimates at low 
temperature and smoothly connects to it at higher tem- 
perature. Thus the range of length scales where the sys- 
tem behaves effectively as a Bragg glass is wider than 
previously expected. We also computed the bare core en- 
ergy dependence of which exhibits stretched exponen- 
tial behaviour. It would be interesting to further explore 
numerically these systems, particularly at weak disorder. 



(13) 



In the limit of large E® this gives: 
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a e ( ^ )3 (14) 

with E* w 2 V /2C(TU. It is noteworthy that a very 
recent numerical workO finds a similar law with an ex- 
ponent ~ 0.7. The above result can be written E® ~ 
(ln^£>) 3 / 2 in agreement with the general scaling of opti- 
mal disorder energies in this problem. 
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